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Magnetic texturing on the surface of a topological insulator allows the design of wave guide 
networks and beam splitters for domain- wall Dirac fermions. Guided by simple analytic arguments 
we model a Dirac fermion interferometer consisting of two parallel pathways, whereby a newly 
developed staggered-grid leap-frog discretization scheme in 2+1 dimensions with absorbing boundary 
conditions is employed. The net transmission can be tuned between constructive to destructive 
interference, either by variation of the magnetization (path length) or an applied bias (wave length). 
Based on this principle, a Dirac fermion transistor is proposed. Extensions to more general networks 
are discussed. 
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Remarkable progress has been made in the understand- 
ing and realization of metallic surface states on three- 
dimensional (3D) topological insulators (TI) pQ. Both a 
theoretical foundation for the characterization of Dirac 
fermion surface states in form of topological invariants, 
as well as subsequent experimental verification of Dirac 
cones of helical surface state have been presented [2HZ]- 
In essence, the ingredients for such TI states are an insu- 
lating bulk material with strong spin-orbit coupling pro- 
ducing a band inversion ( "negative gap" ) within the bulk 
and the conservation of time-reversal (TR) symmetry [Ti- 
ll]. Following a bulk-boundary correspondence principle 
or, alternatively, treating the surface as a domain wall 
where the effective mass changes its sign, one is lead to 
topologically protected gapless surface states of helical 
nature. They can be described by an effective 2D field 
theory which leads to the effective particle Hamiltonian 

H = v (p x a x + Py(Ty) + M(x, y)a z + V(x, y)t (1) 

near the TR degeneracy point and a 2D (2+1) Dirac 
equation for the surface states [4j|8]. Here, v, Pi, and 
denote, respectively, the effective velocity, momentum, 
and Pauli matrices for i = x, y. We have already added a 
space-dependent mass M(x, y) term which is introduced 
by an interaction which breaks TR symmetry, such as the 
proximity of ferromagnetic layers [4] or magnetic doping 
[9]. The order-of-magnitude for the mass-gap is several 
10's of meV. In contrast to graphene, the Pauli matri- 
ces characterize the real spin whereby the spin direction 
S cx z x cr, with z being the surface normal vector, is 
locked perpendicular to the wave vector. This property, 
as well as the presence of a single Dirac cone, on one side 
of the insulator, are promising differences compared to 
graphene when device applications are considered. Up 
to recently, transport experiments have been hampered 
by impurity bands and by problems with the position- 
ing of the Fermi energy inside of the gap of the insulator 
[3 [lOl [TT] . However, more recently this problem seems to 
have been overcome [lOl [TT] . To the best of our knowl- 
edge, experimental results for the coherence length on 



bulk TIs are not yet available, however, measurement of 
Aharonov-Bohm oscillations in the magnetoresistance of 
layered Bi 2 Se3 nano-ribbons have lead to an estimated 
coherence length of around 0.5 fim at low temperatures, 
which justifies our entirely coherent treatment below [T2] . 

Realization of chiral surface states opens new prospects 
for device designs which utilize their coherence proper- 
ties, on one hand, and requires the development of the- 
oretical tools to simulate the dynamics associated with 
them, on the other hand. This Letter addresses both 
issues: first we show by simple analytic argument that 
magnetic structuring on the surface of a topological in- 
sulator allows the formation of chiral channels and beam 
splitters in analogy to their optical counterparts. We 
then propose and model a two-armed quantum interfer- 
ometer, which can be controlled by magnetic structur- 
ing and an electric gate bias. Based on this interfer- 
ometer, we propose a Dirac fermion transistor. Lastly, 
we present an efficient numerical algorithm to evaluate 
the fully time-dependent (ballistic) dynamics of Dirac 
fermion wave packets in 2+1 dimensions, which allows 
a simple inclusion of absorbing boundary conditions at 
the simulation boundaries. This method is used to nu- 
merically demonstrate the time-dependent fermion dy- 
namics within the interferometer under variation of gate 
bias and other structural properties. Generalization to 
more complex networks is discussed. 

It is well established that a ID edge channel forms 
along the domain boundary between two 2D TI surface 
areas when the effective mass changes its sign, see Fig{l] 
[U HI [13]. This edge state is non-degenerate and chi- 
ral in the sense that its spin is locked perpendicular to 
its direction of momentum. X-shaped interconnects be- 
tween edge channels, as sketched in Fig(2j are produced 
by magnetic texturing and provide a "beam splitter" for 
chiral fermions. To understand its principle we recall 
the low energy eigenstate, localized at the effective-mass 
domain- wall. Let, as in Fig(TJ the mass be allowed to 
vary in ^/-direction only, with y' = —xsinO + ycosO, 
and let us assume translational invariance in x'-direction 
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FIG. 1: (color online). Mass domain- wall on the 2D surface 
of a 3D TI: For \im y '^± 00 m(x , ,y') — m± and ra_ra+ < a 
chiral domain wall state localized in ^/-direction forms. 
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FIG. 2: (color online). Wave packet propagation in an asym- 
metric beam splitter consisting of two linear domain walls 
intersecting at 45 degrees with constant masses m+ > and 
m_ (= -m+ ) in region "+"and "— ". The color (brightness 
variation) encodes the phase of the upper spinor component. 
Arrows indicate the allowed direction of propagation. 



(x f = xcosQ + 2/sin0), such that \im y '^± OQ m(y') = m± 
and m_m + < 0. Then a domain- wall channel state 
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(2) 

exists with E = ±vh k' x , with the upper sign for m + < 
and ra_ > 0, and the lower sign for m_ < and ra + > 0. 
Note that a constant potential V in Eq{l] simply adds to 
E. 

Now consider the symmetric junction A in Figj3^a), 
where the domain- wall channel in -f ^-direction splits into 
the two channels along the indirection (9 = ±7r/2). 
At this junction we can express the incoming wave func- 
tion \X, +) as a superposition of the F-out-channel states 
= l/v^(|y,+>+i|-y,+»- The two compo- 
nents |y, +) and |— y,+), respectively, channeled along 



path 1 and 2, necessarily recombine at junction B in 
Figj3ja) to \Xj +) as long as the the phase difference be- 
tween the upper and the lower channel is zero. How- 
ever, the relative phase <j) between path 1 and 2 can 
be set to a desired value by adjusting the gate bias 
V as shown in Figj3^b) and/or by changing the rela- 
tive path length by selectively changing the size of in- 
dividual magnetic domains. Using the decomposition 
|±Y,+) = l/\/2(|±X,+) - i |=F^,+)) one has at junc- 
tion 5, up to an overall phase, 

|V#)> = cos(0/2) |X, +) + sin(0/2) \-X, +) . 

This gives \X,+) for (j) = and for <\> = tt 

which corresponds to fermions propagating, respectively, 
from junction B out to the right (towards drain 2) and 
fermions going to the left (towards junction A). The 
phase difference determines the out-channel. A phase 
shift can be introduced conveniently by an (piece-wise 
constant) electric potential V, which locally shifts the 
dispersion by V, leading to a phase difference <p = 
— V • L/(vh), with L being the length of the path where 
the potential is applied. The transmission probability 
into drain 2 in Figj3|b) may be written 



T 2 = |(X,+|^)| 2 =cos 2 



VL \ 



and the transmission probability into drain 1 in Fig(3jb) 



is 



T 1 = \(-X,+\i>)\ 2 =^ > (^j=l-T 2 . 

Switching gate potentials, rather than the relative size 
and shape of magnetic domains, most likely is the more 
convenient way to control the device. Gate potentials can 
also be used to compensate imperfect magnetic texturing 
in experimental realizations. 

The above pedagogical discussion was limited to (the 
asymptotic regions of) rectangular symmetric junctions 
of domain- wall states, which correspond to symmetric 
beam splitters in optics, albeit, without back scattering 
channel. Both the selection of asymmetric mass profiles 
and junctions where the two in- and two out-channels 
meet at arbitrary angles can be used to design asymmet- 
ric beam splitters as shown in Figj2j These are best mod- 
eled numerically. Studies of Dirac fermions traditionally 
have been conducted by the high-energy physics com- 
munity, however, recently both atomic and condensed 
matter theorists have picked up on this topic as well 
P~4j [15]. Fermion doubling has been a longstanding 
problem in lattice simulations using local discretization 
schemes [16] . Here we propose a staggered-grid leap-frog 
discretization scheme shown in Fig(4] and use the nota- 
tion rj)(xj,y k ,t n ) = tl>j ik = (uV k ,v% k ), where n and j,k, 
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FIG. 3: (color online). Coherent Dirac fermion transistor 
on a TI surface: Chiral channel states form at the domain 
boundaries between magnetization regions of opposite direc- 
tion (M| /M\). There are two allowed paths trough the 
structure (wavy lines). When a gate voltage is applied the 
paths pick up an additional phase leading to destructive or 
constructive interference depending on its magnitude and the 
device can be switched from transmission to the left (drain 1) 
to transmission to the right (drain 2). 



respectively, are the discrete time and space indices: 
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A stability and dispersion analysis for this linear system 
for m, V G 5R, by Fourier analysis, shows that the norm 
\\1pW2 is preserved. For At = Ax this method has the 
dispersion shown in Figj5} providing the exact relativis- 
tic dispersion along the x- and y-axis, while one extra 
Dirac cone is shared by the 4 corners of the Brillouin 
zone. At the expense of perfect dispersion along the x- 
and y-axis and breaking chiral-symmetry one can get a 
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FIG. 4: (color online). Leap-frog staggered-grid scheme: The 
left part shows the time- stepping where 1) the new u compo- 
nents (blue) are computed by the previous u and the spatial 
differences of old v- values. 2) Then (knowing u at t n +i) the 
new v at t n +i (red) are computed. The mass m and potential 
V enter the scheme in a Crank-Nicolson-type time averaging 
over the current- and previous-time values. The right part of 
this figure shows the scheme for the spatial derivatives. 




FIG. 5: (color online). Dispersion relation for the leap-frog 
staggered-grid scheme in normalized units and At = Ax = 
Ay = a, m — 0: The contour lines encode normalized energy 
e • a. It shows perfect linear dispersion in x- and y- direc- 
tion thus avoiding phase and group velocity errors, essential 
for observing the correct interference behavior for rectangular 
domain- wall junctions. 



single Dirac cone by adding a Wilson mass term lead- 
ing to staggered Wilson fermions [17 j. Our numerical 
scheme also allows an easy implementation of absorbing 
boundary conditions by putting a layer with imaginary 
potential around the simulation region which avoids spu- 
rious reflections at the simulation boundaries or the need 
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FIG. 6: (color online). Wave propagation in the interferom- 
eter for (a) phase difference = and (b) phase difference 
(f) — 7r. The wave packet is shown for increasing time as it 
propagates along the domain wall channels. The color (or 
brightness variation) encodes the probability density \ijj\ 2 ' . 
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FIG. 7: (color online). Transmission of the coherent fermion 
transistor (to drain 2) as a function of the gate voltage: (a) 
time resolved transmission for different gate voltages; (b) 
transmission as a function of gate voltage. The solid line 
represents the analytic expression, and the circles are results 
of the simulation. 



wave packet is prepared with mean energy of 30 meV. 
The mass gap due to the magnetization is chosen to be 
100 meV. For a phase difference cj> = in part (a) of FigjH] 
the wave packet leaves the interferometer to the right, 
whereas for <ft = 7r in part (b) the wave packet is forced 
to the left back into the interferometer. Note also that 
the wave packet stays compact due to a faithful represen- 
tation of the Dirac fermion dispersion in our simulation. 
The time resolved transmission for the simulated wave 
packet and varying gate bias is shown in Fig{7^a). Over 
the range of 30 meV the out-channel to drain 2 can be 
tuned from open to closed. The interference fringes ver- 
sus gate bias for transmission probability T2, plotted in 
FigjT^b), confirms the analytic prediction. 

Finally, let us conclude by summarizing the main 
points of this Letter and give an outlook for further re- 
search. We propose electronic networks of chiral domain- 
wall states, arising from magnetic domains built into the 
surface of topological insulators, in analogy to optical 
wave guides. We propose and analyze numerically a 
Dirac fermion beam splitter. Combination of two such 
splitters leads to an interferometer which can be con- 
trolled by an electric gate and extended to produce a 
"Dirac fermion transistor" . The extension to more com- 
plex Dirac fermion networks is straight-forward in princi- 
ple, whereby, curved domain walls serve as one-way wave 
guides and intersections of domain wall boundaries pro- 
vide beam splitters, phase differences can be controlled 
by electric gates, and charge can be provided by electric 
contacts. 

We acknowledge support from the Austrian Science 
Foundation project I395-N16. 



for periodic boundary conditions. 

Figj2] shows results for the wave packet propagation 
through an asymmetric beam splitter made up of two 
linear domain- walls intersecting at 45 degrees. Areas 
marked by "+" and "— ", respectively, correspond to re- 
gions with constant mass m + > and m_= -m + . Unlike 
for rectangular junctions, the splitting ratio also depends 
on the shape of the incident wave packet. This can be un- 
derstood by accounting for the transverse spatial extent 
which, see Eq(2j is controlled by the mass profile of the 
wave function associated with the channel states lead- 
ing to their mutual overlap near the junction. Therefore, 
both, mass profiles and the angle of intersection influence 
the properties of the beam splitter. A more detailed anal- 
ysis of asymmetric junctions will be given elsewhere. 

Results for the interferometer, laid out in Figj3J are 
given in Figs (6] and [7| A gate bias is applied on the up- 
per half-plane (path 1) over a channel length of ~ 260 
nm. The distance source-drain is ~ 400 nm. The initial 
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